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We study some properties of the Ising model in the plane of the complex (energy/temperature)- 
dependent variable u — e~ 4K , where K = J/(fcsT), for nonzero external magnetic field, H. Exact 
results are given for the phase diagram in the u plane for the model in one dimension and on infinite- 
length quasi-one-dimensional strips. In the case of real h — H/(kBT), these results provide new 
insights into features of our earlier study of this case. We also consider complex h — ///(fcsT) and 
/i = e~ 2h . Calculations of complex-M zeros of the partition function on sections of the square lattice 
are presented. For the case of imaginary h, i.e., = e l ° , we use exact results for the quasi- ID strips 
together with these partition function zeros for the model in 2D to infer some properties of the 
resultant phase diagram in the u plane. We find that in this case, the phase boundary B u contains 
a real line segment extending through part of the physical ferromagnetic interval < u < 1, with 
a right-hand endpoint u r he at the temperature for which the Yang-Lee edge singularity occurs at 
(i — e ±l ° . Conformal field theory arguments are used to relate the singularities at u r h e and the 
Yang-Lee edge. 
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I. INTRODUCTION 

The Ising model serves as a prototype of a statistical 
mechanical system which undergoes a phase transition 
in the Z2 universality class with associated spontaneous 
symmetry breaking and long-range order. At tempera- 
ture T on a lattice A in an external field H , this model is 
defined by the partition function Z = } e ~^ H ^ with 
Hamiltonian 



H 
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where <jj — ±1 are the classical spin variables on each 
site j S A, (3 = (fc^T) -1 , J is the spin-spin exchange 
constant, and (jf) denote nearest-neighbour sites. We 
use the notation K — (3J, h = f3H, 



-4K 



-2h 



(1.2) 



The free energy is F = —ksTf, where the reduced free 
energy is / = linin^oo n^ 1 InZ, with n being the num- 
ber of lattice sites in A. Physical realizations of the 
Ising model include uniaxial magnetic materials, struc- 
tural transitions in binary alloys such as (3 brass, and the 
lattice-gas model of liquid-gas phase transitions. The 
two-dimensional version of the model was important 
partly because it was amenable to exact solution in zero 
external magnetic field and the critical point was charac- 
terized by exponents that differed from mean-field the- 
ory (Landau-Ginzburg) values Q-p^- Although we shall 



phrase our discussion in the language of the Ising model 
as a magnetic system, the results have analogues in the 
application as a lattice-gas model of a liquid-gas phase 
transition, in which e~ 2h corresponds to the fugacity. 

Just as one gains a deeper understanding of functions 
of a real variable in mathematics by studying their gen- 
eralizations to functions of a complex variable, so also it 
has been useful to study the generalization of h and K 
from real to complex values, as was pioneered by Yang 
and Lee 0] and Fisher respectively. In this con- 
text, one finds that the values of u where the model has 
a paramagnetic-to-ferromagnetic (PM-FM) phase tran- 
sition and a paramagnetic-to-antiferromagnetic (PM- 
AFM) phase transition occur where certain curves in the 
complex-u plane cross the positive real-u axis. These 
curves define boundaries B u of complex-it extensions of 
the physical phases of the model and arise via the accu- 
mulation of zeros of the partition function in the thermo- 
dynamic limit. One of the aspects of complex- it singular- 
ities studied in early work was their effect on the conver- 
gence of low-temperature series expansions [6] . Although 
the two-dimensional Ising model has never been solved 
exactly in an arbitrary nonzero external magnetic field 
H , the free energy and magnetization have been calcu- 
lated for the particular imaginary values h = i(2£+l)7r/2 
with £ £ Z, which map to the single value jj, = — 1 
0) S I n previous work we presented exact deter- 
minations of the boundaries B u for this fj, = — 1 case on 
the square, triangular, and honeycomb lattices, as well 
as certain heteropolygonal lattices d, [ljl H3 ■ We inves- 



tigated the complex-it phase diagram of the Ising model 
on the square lattice for physical external magnetic field 
in (lol | , using calculations of partition function zeros and 
analyses of low-temperature, high- field (small- |it|, small- 
series to study certain singularities at endpoints of 
lines or curves of zeros. In that work we considered real 
h and the complex set h = h r + Htt/2 (where t € Z) 
that yield real \i. (Our notation [i — e~ 2h follows that in 
the original papers on series expansions [H, that we 
used in [l(| and should not be confused with the chemical 
potential in a liquid-gas context.) 

In this paper we continue the study of the complex-w 
phase diagram of the Ising model for nonzero h. We 
present exact results for lattice strips, including their 
infinite-length limits, and calculations of partition func- 
tion zeros on finite sections of the square lattice. For 
the special case of real h and the subset of complex h 
that yield real pi, our exact results for these strips pro- 
vide new insight into properties that we found for the 
2D Ising model in [l0( . Complex-w zeros of the 2D Ising 
model partition function with nonzero field have been 
studied further in subsequent works [3, EH- Among 
general complex values of h, we pay particular attention 
to the case where h is pure imaginary. This is of inter- 
est partly because of an important property of the Ising 
model that was proved by Yang and Lee [3| , namely that 
for the ferromagnetic case (J > 0), the zeros of the parti- 
tion function in the y, plane lie on the unit circle \fi\ = 1, 
i.e., correspond to imaginary h. In the limit where the 
number of sites n — ► oo, these zeros merge to form the 
locus B^ comprised of a connected circular arc /i = e 10 , 
where id = —2h, passing through fi = — 1 (i.e., 9 = tt) 
and extending over on the right to a complex-conjugate 
pair of endpoints at e ±t6c . This result applies for the 
Ising model in any dimension; indeed, it does not require 
A to be a regular lattice. One interesting question that 
we address is the following: what is the phase boundary 
B u in the u plane for /i = e 10 when 9 is not equal to 
one of the two exactly solved cases, i.e., 9 =^ mod w. 
We answer this question with exact results for quasi-lD 
strips and study it with partition function zeros for 2D. 
We find that, in general, B u contains a real line segment 
extending through part of the physical ferromagnetic in- 
terval < u < 1, with a right-hand endpoint u r h e at 
the temperature for which the Yang-Lee edge singularity 
occurs at /i = e ±t6 . We use conformal field theory argu- 
ments to relate endpoint singularities in the u plane for 
real and imaginary h to the Yang-Lee endpoint (edge) 
singularity. 



II. RELEVANT SYMMETRIES 

We record here some basic symmetries which will be 
used in our work. On a lattice with even (odd) coordi- 
nation number, the Ising model partition function Z is 
a Laurent polynomial, with both positive and negative 
powers, in u (in ^/u). Z is also a Laurent polynomial in 



fi and hence, without loss of generality, we consider only 
the range 

--< lm(h) < - . (2.1) 

Furthermore, Z is invariant under the simultaneous 
transformations h — > —h, Uj — ► — <7jVj G A. The sign 
flip h —>■ — h is equivalent to the inversion map 

1 

h -> -h <-> //—> — . (2.2) 
A 4 

Hence, in considering nonzero real h, one may, with no 
loss of generality, restrict to h > 0. More generally, in 
considering complex h, one may, with no loss of general- 
ity, restrict to the unit disk in the /i plane, \fj,\ < 1. It is of 
particular interest to consider two routes in the complex 
fi plane that connect the two values of /i where the 2D 
Ising model has been exactly solved, viz., /i = 1 (h = 0) 
and n = — 1 (h = iir/2). The first such route is the one 
that we used in viz., the real segment — 1 < fj, < 1. 
A second route proceeds along the unit circle [i — e %e . In 
view of the above symmetries, it will suffice to consider 
this route as 9 increases from to 7r. If /i £ R, then 
the set of zeros of Z in the u plane is invariant under 
u — > u* and hence the asymptotic locus B u is invariant 
under u — > u* . 

The invariance of the set of complex-temperature zeros 
in u under the complex conjugation u — > u* holds not 
just for real \x but more generally for /x on the unit circle 
= 1. This is proved as follows. For any lattice A, 
Z(A;u,li) = Z(A;u,l/n). Now if and only if fi = e i0 
(with real 9), then /i -1 = fi*. Hence, for this case of 

f j, = e ie Z(A;u,n) = Z(A;u,fi*) . (2.3) 

Now since Z(A, u, n*) = [Z{ A, u*, //)]*, it follows that if 
^ = e ie (with 9 G K), then 

Z(A;u,e M ) = ^ Z(A; u*, e tf> ) = , (2.4) 

so that the set of zeros of the partition function in the 
u plane is invariant under complex conjugation for this 
case. 

III. PROPERTIES OF ID SOLUTION 
A. General 

Because of its simplicity and exact solvability, the 
ID Ising model provides quite useful insights into prop- 
erties for complex temperature and field. As is well 
known, the Perron-Frobenius theorem guarantees that 
the free energy and thermodynamic quantities of a spin 
model with short-ranged interactions are analytic func- 
tions for any finite temperature on infinite-length strips 
of bounded width. The physical thermodynamic prop- 
erties of the Ising model on quasi-one-dimensional strips 



2 



are thus qualitatively different from those on lattices of 
dimensionality d > 2. However, in addition to the physi- 
cal critical point at zero temperature, switching of dom- 
inant eigenvalues of the transfer matrix and associated 
non-analyticity in these quantities can occur for complex 
K and/or h. Indeed, the properties of the model, on 
quasi-lD strips, for complex u and [i exhibit some inter- 
esting similarities to those on higher-dimensional lattices. 
For example, the phase boundary for the zero-field Ising 
model on the square, triangular, and honeycomb lattices 
exhibits a multiple point at u = — 1. (Here, the term 
"multiple point" is used in the technical sense of alge- 
braic geometry and is defined as a point where two or 
more branches of the curves comprising this boundary 
cross each other.) This feature of the model on the 2D 
lattices is also found to occur for quasi- ID strips such as 
the L y = 2 strip s of the square [l6| , triangular [l7| , and 
honeycomb [18j lattices. Furthermore, the fact that the 
circle theorem of 0] applies in any dimension means that 
there is particular interest in using quasi- ID strips to ob- 
tain exact results on the singular locus B u corresponding 
to a point on the unit circle in the [i plane. 



B. Calculation of B u and Analysis of 
Thermodynamic Quantities 

We begin our analysis of quasi-one-dimcnsional lattice 
strips with the ID line with periodic boundary condi- 
tions, i.e., the circuit graph, C n . The well-known transfer 
matrix is 



e K+h £ -K 



-,-K 



The eigenvalues of this transfer matrix are 



K 



sh{h) ± (sinh 2 (/i) 



+ e" 



-4K 



1/2 



(3.1) 



, (3.2) 



where the + and — signs apply for j = 1 and j ' = 2. The 
eigenvalues have branch-point singularities at u = u e , 
where 



sinh 2 (/i) 



(u + ^-2) 



(3.3) 



(the subscript e denotes "endpoint"). Note that 
det(Ti£>) = 2sinh(2A"), independent of h. The partition 

function is Z(C n ) = Tr[(T 1D ) n ] = £j = i(Air>,j) n - We re- 
strict to even n to avoid frustration in the antiferromag- 
netic case. The reduced free energy is / = ln(\iD,m,axX 
where \\D,dom. denotes the maximal eigenvalue. Equiv- 
alcntly, 



1 



/ = ~ / 7T ln (Aid.i) 2 + (Aid, 2 ) 2 -2A 1Z5 . 1 A 1Di2 cos((/)) 
2 Jo 2ir l 

(3.4) 



i.e., explicitly, 



hi 



cosh(2h)+u+(u— 1) cos(<^>) 



(3.5) 

The complex-M phase boundary B u is the locus of so- 
lutions in u to the condition that the argument of the 
logarithm in eq. (|3.5p vanishes, i.e., 



B u : cosh(2/i) + u + (u - 1) cos(^) = 



(3.6) 



In accordance with the general discussion of symme- 
tries given above, B u is symmetric under h — > —h, i.e., 
(i — > If and only if /i — ±1, it is also symmetric un- 
der K — > —K. i.e., u — > 1/u. The condition (|3.6[) is con- 
dition that there is degeneracy in magnitude among the 
dominant (and here, the only) eigenvalues of the transfer 
matrix, IAlDji] = IAld^I- For real h, this is equivalent 
to the condition that the argument of the square root in 
eq. (|3.2|) is negative. The locus B u is thus a semi-infinite 
line segment on the negative real axis, 



B u 



u < u e 



for real h 



(3.7) 



The right-hand endpoint of this line segment, u r h e = u e , 
occurs at 14 = if and only if /i = 1, i.e., h = 0. As 
\h\ increases, u e moves to the left along the negative real 
axis. For h = 0, the phase boundary B u is noncompact in 
both the 14 and 1/u planes (one implying the other by the 
invariance of B u under the inversion map u — » 1/u), but 
for h 7^ 0, it is noncompact in the u plane but compact 
in the 1/u plane. 

For complex h = h r + in/2, or equivalcntly, —1 < fi < 
0, the term cosh(/i) in the eigenvalues (|3.2[) is imaginary, 
so the condition that these eigenvalues be equal in mag- 
nitude is the condition that the square root should be 
real. Hence, 



B u 



u > u P = 



2 + |mI + ImI" 



for fi < 



(3.8) 



This is a semi-infinite line segment on the positive real 
axis in the u plane with left-hand endpoint uihe = u e . 
For this case of negative real jji, uihe > 1, and uihe —> 1 + 
as fi + 1 — ► + . As \h\ — > oo, upj le — * oo. Again, in 
the 1/u plane, this is a finite line segment from to the 
inverse of the right-hand side of eq. (|3 . 8[) . 

We next determine B u for /i = e 10 on the unit circle 
in the [i plane. In this case, with it e = — sinh 2 (/i) = 
sin 2 {0/2), one has 



B u : u< sin 2 (6/2) for fj, = e l 



(3.9) 



This is a semi-infinite line segment whose right-hand end- 
point u r hs — Ue occurs in the physical ferromagnetic in- 
terval < u r he 5: 1; increasing from u = at 9 = to 
U = 1 as 6 approaches ir from below. For all values of \x 
' on the unit circle except for the points /it = ±1, the locus 
B u is not invariant under u — > 1/u. 
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Finally, for fi = — 1, the eigenvalues are equal in mag- 
nitude and opposite in sign so that, with n = 2£ even, 
Z = 2z (1 — it) . (If n were odd, then Z would vanish). 
Since we keep n even, Z vanishes only at the point u = 1, 
and B u degenerates from a one-dimensional locus to the 
zero-dimensional locus at u = 1. 



C. Singularities at u e 

The physical singularities of the zero-field ID Ising 
model at T = are well known; taking J > without 
loss of generality, the model exhibits exponential diver- 
gences x — /3e 2K and £ ~ (l/2)e 2K in the susceptibility 
and correlation length, together with a jump discontinu- 
ity in the spontaneous magnetization. Here we focus on 
the singularities in thermodynamic quantities at u e for 
nonzero ft. The internal energy per site is 



U 



-J 



1- 



2e 



-4K 



cosh(/i) + yj sinh 2 (ft) + u ^ ^/ sinh 2 (ft) + u. 

(3.10) 

As is evident from this or from an explicit calculation of 
the specific heat Ch = (dU/dT) H , for nonzero ft, h ^ 
i(2£ + 1)tt/2 with £ £ Z, C# diverges at u — u e with 
exponent a' e = 3/2. (We use primes on a and 7 for 
phases with either explicitly or spontaneously broken Z 2 
symmetry.) Applying the standard scaling relation 2 — 
a' = d/yt (where y t is the thermal exponent) at it = u e 



yields yt = 2 at this singularity. For ft 
i.e., fi = — 1, eq. (|3.10p reduces to 



-1) 



J 



so that 



C*(/i=-l) 



tanh(2X) 



2k B K 2 
sinh 2 (2if) 



i(2^ + 1)tt/2, 



(3.11) 



(3.12) 



Hence, if [i = —1, whence u e — 1, i.e., -RT e = 0, the 
specific heat is finite at u e , and a' e — 0. This is the same 
value that we found for the 2D Ising model at fj, = —1, 
U s — 1 in Q, as discussed further below. 
The per-site magnetization is 



M = 



sinh(/i) 



\J sinh 2 (/i) + u 



(3.13) 



For h ^ this diverges at u — u e with exponent (3 e 
— 1/2. The susceptibility per site, x = dM/dH, is 



(3u cosh(/i) 
(sinh 2 (/ l )+ U ) 3 / 2 



(3.14) 



For h 7^ and ft, 7^ i(2^ + l)7r/2, this diverges at it = u e 
with exponent 7^ = 3/2. Applying the scaling relation 
/3 + 7' = 2/ft/j/i (where yh is the magnetic exponent) at 



u e and substituting y t = 2 then yields = 2 at this sin- 
gularity, so that y t — yh at u e . For h = i(2£+ l)n/2, the 
cosh(ft) factor causes x to vanish identically, so that no 
exponent 7^ is defined. Thus, this exactly solved model 
shows that, just as the value ft, = is obviously special 
since it preserves the Z2 symmetry, so also the values 
h = i(2£ + l)n/2 are special, leading to different values 
of singular exponents at u e than the values at generic 
nonzero values of ft. 

We denote the density of zeros on B u as g(u). As 
u approaches a singular point u s on B u , the density of 
zeros is related to the critical exponent for the specific 
heat 0,IH 



9{u) 



U - M, 







(3.15) 



(This exponent would be denoted 1 — a if it applies at 
the critical point, as approached from the physical high- 
temperature phase.) The singular point u s may be the 
critical point, u Cl as in the case h — 0, or the arc end- 
points u e and u* studied in [To| in the presence of a real 
nonzero field, and, in our present discussion we are inter- 
ested in the endpoint of the line segment at u e . 

Let us first consider the boundary B u for physical ft. 
This locus is the solution to eq. (|3.6p and the density 
of zeros is proportional to d<p/du. To begin, we consider 
real ft. It is convenient to introduce a positive variable 
u! = — u. If one normalizes the density according to 



then 



9i u ') = 



du' g(u') = 1 , 



2 cosh(ft) 



tt(1 + u')Ju' - sinh 2 (ft) 



(3.16) 



(3.17) 



In the neighborhood of a point where the free energy is 
singular, one can write, as was done in [lol ]. 



(3.18) 



From the discussion above, one already knows for ft 7^ 
0, a' e — 3/2 at u e except if ft = (21 + l)wr/2, where 
a' e = and g(u) vanishes identically, reflecting the above- 
mentioned fact that B u degenerates to a point at u — 
u e = 1 . These findings are in agreement with the present 
analysis of the density of zeros; for ft 7^ (21 + l)wr/2, 
expanding eq. (|3.17p as vl — u' e — > + , we have g(u') — > 
{2/ir)/y/u'-u' e , so 1 - a' e = -1/2, i.e., a' e = 3/2. 

We next show the close relation between this singular 
behavior of the density of zeros on B u as one approaches 
the endpoint u e with the singular behavior of the zeros 
on B^ as one approaches the endpoint [i e of that locus. 
We focus on the case J > and ft imaginary, for which 
By is an arc of the unit circle fi — e %e extending clockwise 
from 9 = 7r to 9 = 9 e and counterclockwise from 9 = tt 



4 



to 9 = —9 e . The density of zeros on i3 M , denoted g(0), 
has the singular behavior at the endpoint e 10 " given by 

g{9) - (6> - 6 e ) a as 9 - 9 e -► 0+ . (3.19) 

With i0 = -2h, eq. ([376]) becomes 

cos0 + u+(w-l)cos^ = O . (3.20) 

Letting <^> range from to 2tt, one sees that the endpoints 
occur at 

8 e = arccos(l — 2u) = 2arcsin(y / u ) , (3-21) 
i.e., in terms of /j,, 

// e ,/i* = 1 - 2u ± 2iv/u(l - u) . (3.22) 



The density of zeros, i.e., the number of zeros N z between 
9 and 6* + d9, is given by differentiating with respect to cf> 
and noting that the totality of these zeros corresponds to 
the range < d>< 2n: g(9) = dN z /d9 = {2n)- 1 c 
The density is Q 



g{0) 



sin(0/2) 



(3.23) 



2%Jsm 2 {9/2) - u 



for sin 2 (6>/2) > u and g{9) = for sin 2 ((9/2) < u. 

This density diverges as 9 — 9 e -> 0+, with the Yang- 
Lee edge exponent a = —1/2 0, [32J. Given the scaling 
relations 



d-2 



d + 2- 



H = — - i 
»7 Vh 



(3.24) 



the result cr e = —1/2 is equivalent to yh = 2 at u e . As 
was noted in for the antiferromagnet (J < 0), the 
zeros in the /i plane form a line segment on the nega- 
tive real [i axis. The singularities in the density g{p) at 
the endpoints of this line segment are again square root 
singularities. Thus, for this exactly solved ID model, 



1 



a, = a = 



-- (ID) 
2 y ' 



(3.25) 



That is, the exponent 1 — a' e describing the singular be- 
havior in the density of partition function zeros in the 
locus B u in the u plane as one approaches the endpoint 
u e of this locus is the same as the exponent a = —1/2 
describing the singular behavior in the density of zeros 
in the locus as one approaches the endpoints of this 
locus in the \i plane. This shows, as we have emphasized 
in our earlier work [ToL fll| , the value of analyzing the 
singular locus £>, including its slice B u in the u plane for 
fixed fj, and its slice £> M in the /i plane for fixed u, in a 
unified manner. Indeed, the value of such a unified ap- 
proach to this singular locus was recognized in general in 
early works such as [2l|, [22j ■ We note also that for the 
ID Ising (and Potts) models, there is a duality relation 
connecting temperature and field variables [23, 24]. The 



intertwined relation of the two relevant variables K and 
h for a (bi)critical point is at the heart of the analysis of 
the scaling limit (T — T c ) — * 0, H — > in terms of the 
scaling variable h/(T — T c ) yh ^ Vt , where yt and yh denote 
the thermal and magnetic exponents. A difference is that 
in the present analysis, the singular point (s) u s is (are) 
not, in general, the physical critical point. 

For the case of complex [i with ^ 1, our analysis 
of the singularity at u e goes through as before. However, 
in this case, because it is no longer true that yT 1 = 
fi*, the coefficients of the powers of u in the Laurent 
polynomial comprising Z are not real, so the set of zeros 
in the u plane for a given /u, is not invariant under complex 
conjugation. 



IV. EXACT SOLUTION FOR TOROIDAL 
LADDER STRIP 

A. General Calculation 

Here we consider the ladder strip of the square lattice 
with doubly periodic (i.e., toroidal) boundary conditions. 
These boundary conditions have the advantage of mini- 
mizing finite-size effects. They also have the merit that 
for any length L x , all of the sites on the lattice have the 
same coordination number, equal to the value of 4 for the 
infinite square lattice. The periodic transverse boundary 
conditions entail a double bond between the sites on the 
upper and lower sides of the ladder. For a given length 
L Xl the strip has n = 2L X sites. In the basis 



{(:)• 



:)} 



(4.1) 



the transfer matrix for this toroidal ladder (t£) strip is 



The determinant is 



/ 


AK+2h 


e h e h 1 




' e h 


1 e~ iK e~ h 




e h 


e -AK x e -h 


V 


1 


e -h e -h e AK-2h 



(4.2) 



det(Ttf) = [2sinh(2if)] 



4_ (l-«) 4 



(4.3) 



evidently independent of h. The partition function is 

4 

^ = Tr[(r«)^] = 5>^) n/2 (4.4) 
i=i 

where 

A«,i =l-u (4.5) 
and the three other \tt,j 's are roots of the cubic equation 
A 3 + a 2 A 2 + aiA + a = (4.6) 
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where 



a 2 = -(l + u + u 1 {n + n 1 )j 



a i 



(1-u) 1 + u{fi + fi- 1 + 1) 



and 



a = 



(u - If 



(4.7) 



(4.8) 



(4.9) 



The reduced free energy is given by / = (1/2) ln(A«, mQX ). 
Because of the cumbersome form of the solutions to the 
cubic equation, we do not display the explicit results for 
thermodynamic quantities such as the specific heat, mag- 
netization, and susceptibility. Our primary purpose in 
the analysis of this toroidal strip is to determine the 
boundary B u for a given fi and to glean some insights 
from exact results on this boundary for the case of the 
model in 2D. 
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FIG. 1: Complex- temperature phase boundary B u and partition 
function zeros in the u plane, for the Ising model with h = 0, 
i.e., n = 1, on an L y = 2 strip of the square lattice with toroidal 



boundary conditions. Zeros are shown for L x 



200. 



B. Properties at Some Special Points 

In general, a point (u, fi) is contained in the singular 
locus B u if there is a switching of dominant eigenvalues 
of the transfer matrix. We can thus immediately derive 
some results on this locus by considering some special 
cases. 

For u = 1, the eigenvalues are Xuj = for j = 1, 2, 3 
and \ti,4 = (1 + /i) 2 //i, so the eigenvalues are equal at 
this point if and only if /i = — 1. For this value, they all 
vanish, as does the partition function. Hence, the zero of 
the partition function at this point has a multiplicity of 
n/2. For the quasi-lD strips considered here, this point 
(it, /i) = (1, —1) occurs where B u degenerates to a point. 
In contrast, for the square lattice, it is contained as part 
of a one-dimensional locus B u [§]. 

For u = — 1, the eigenvalues are Xu,j = 2, j — 1,2 and 

A«,i = [(1+M) 2 ± V(v - 1) 2 (1 + 6^ + A* 2 ) ] (4-10) 

where the ± sign applies for j — 3,4. For fi — ±1, all 
four of these eigenvalues have magnitudes equal to 2, so 
the points 



-1,±1) S B 



(4.11) 



For real /x > 0, fi ^ 1, \Xa^\ is smaller than 2, de- 
creasing to as jj, — * or /j, — > oo, while \Xte,<k\ is larger 
than 2, approaching infinity as [i — > or /i — > oo, so that 
there are no further switchings of dominant eigenvalues 
for these values of /i. We next consider the real interval 
(i < 0. The polynomial in the square root in eq. ()4.10[) is 
negative for -(3 + 2^2) < fi < -(3-2\/2) and \Xu,j\ = 2 
for all four j = 1, 2, 3, 4 for this interval. Hence, 



B u D {u = -1} for 



(3 + 2V2) < n < -(3 - 2V2) 
(4.12) 



for this strip. Although we give the full range of /z, we 
recall that, owing to the /i <-> symmetry, it is only 
necessary to consider the interior of the disk = 1 since 
the behavior of B u determined by \fj,\ in the exterior of 
this disk is completely determined by the values of /i in 
the interior. 



C. n = 1 

We now proceed with our analysis of the complex-u 
phase diagram for the infinite- length limit of this toroidal 
ladder strip for specific values and ranges of /x. For the 
zero-field case (j, = 1, the three eigenvalues in addition to 
A«,i, are 



At£,2 — 



1 



and 



u + u 1 + 2 ± \fuF 
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(4.13) 



(4.14) 



where the ± sign applies for j = 3, 4, respectively. For 
this case, under the symmetry transformation K — > 
— K, the first two eigenvalues are permuted according to 
A«,i — > -A«, 2 , A t £ j2 — > — A«,i, while the last two, Xu.3 
and A t £,4, are individually invariant. 

For this h = case, in the limit L x — > 00, the bound- 
ary B u consists of an inner closed curve shaped like a 
lima bean passing through the origin u = where it 
has an involution, and through the point u = — 1. The 
rest of B u , which is related to this inner part by the 
u — * 1/u symmetry, passes through u = — 1 and extends 
to u = ±ioo. The point u = —1 is an multiple point 
of osculation type, where the inner and outer curves on 
B u coincide with equal (vertical) tangent. The locus B u 
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FIG. 2: Complex-temperature phase boundary B u and parti- 
tion function zeros in the u plane, for the Ising model with 
h = (1/2) In 2, i.e., fi = 1/2, on a ladder strip with toroidal bound- 
ary conditions. Zeros are shown for L x = 200. 

thus separates the u plane into three regions, which in- 
clude the respective three intervals of the real axis: (i) 
R\\ u > 0, where \u,z is the dominant eigenvalue; (ii) 
i?2: u < — 1, where Xu.i is dominant; and (iii) R3: 
— 1 < u < 0, where \te,2 is dominant. Thus, the outer 
curve is the solution locus of the equation |A«,i| = |At^,3|, 
while the inner bean-shaped curve is the solution locus of 
the equation (A^^l = lA^^j. The outer curves cross the 
imaginary axis at u — ±(\/2 + 1)?, while the inner curves 
cross at the inverses of these points, u = t(v / 2 — 
In Fig. [T] we show a plot of complex-temperature zeros 
calculated for a long finite strip, which clearly indicate 
the asymptotic locus B u . 



D. < n < 1 

We next consider nonzero h, recalling that, owing to 
the symmetries of Z under h —* —h, we can, without 
loss of generality, restrict to the rim and interior of the 
unit disk < 1 in the fi plane. As h increases from 
zero through real values, i.e., /i decreases from 1, the 
part of the locus B u that passed through u — for h = 
breaks apart into two complex-conjugate arcs whose end- 
points move away from the real axis. The outer curves on 
B u continue to extend to infinity in the u plane, passing 
through the origin 1/u = of the e AK plane. This is a 
consequence of the fact that a nonzero (finite) external 
magnetic field does not remove the critical behavior asso- 
ciated with the zero-temperature PM-AFM critical point 
of the Ising model on a bipartite quasi-one-dimensional 
infinite-length strip. The locus B u continues to intersect 
the negative real axis, at the point 

u = . (4.15) 



The outer part of the locus B u is comprised of two 
complex-conjugate curves that extend to complex infin- 
ity, i.e. pass through 1/u = 0. The locus B u separates the 
u plane into two regions: (i) region R\, which contains 
the real interval — /i -1 < u < 00, where the root of the cu- 
bic with greatest magnitude is dominant, and (ii) region 
i?2, which contains the real interval —00 < u < —fi , 
where Xte.i is dominant. The region R3 that was present 
for h = is no longer a separate region, but instead is 
contained in R± . As an illustration of the case of nonzero 
h, we show in Fig. [2] a plot of the phase diagram for 
/i = 1/2, for which B u crosses the real u axis at u = — 2. 
For this value of fi, the arc endpoints on B u are located 
at u ~ 0.149480 ± 0.376522i, which are zeros of the poly- 
nomial 

64u 8 + 128u 7 + 1252u 6 + 1864u 5 + 3448u 4 - 1060u 3 + 
+937u 2 - 108u + 36 , (4.16) 

which occurs in a square root in the solution of the cubic 
equation (|4.6[) . 

The behavior of this exactly solved example provides 
a simple one-dimensional model of the more complicated 
behavior on the square-lattice. For the 2D case with 
any nonzero h, the part of the singular locus B u that 
intersected the real u axis for h = at the position of 
the PM-FM critical point, upm-fm = 3 — 2^2, breaks 
open, with the two complex-conjugate endpoints moving 
away from the real axis, as shown in Fig. 4 of [10(. This 
breaking of the boundary and retraction of the arc end- 
points away from the point upm-fm is in accord with a 
theorem that for nonzero (physical) H, the free energy is 
a real analytic function in an interval from T = beyond 
T c for the PM-FM transition, i.e., in this case, from u = 
along the real u past the point u = u c (25|. For our ex- 
actly solved quasi-lD strips, the PM-FM critical point is 
at u = 0, which is thus the analogue of upm-fm- So the 
motion of the right-hand endpoint of the semi-infinite line 
segment in eq. (13. 7[) . moving left, away from the point 
u = 0, as h increases in magnitude from zero (through 
real values), is analogous to the motion found in [l(| of 
the arc endpoints away from the real axis. The ladder 
strip exhibits a behavior (shown in Fig. [5]) even closer to 
that which we found in the 2D case, namely the breaking 
of the curve on B u that passes through the former crit- 
ical point and the retraction of the complex-conjugate 
endpoints on B u from the real axis. 

E. m = -1 

We can also use our results to consider the complex- 
field value fj, = — 1 and the interval — 1 < fi < 0. We 
begin with the value fj, = — 1. Here the eigenvalues of the 
transfer matrix take the simple form Xte,i = 1 — u as in 
eq. (|4.5[) and, for the three others: 

Xaa^l-u' 1 (4.17) 
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FIG. 3: Complex- temperature phase boundary B u and partition 
function zeros in the u plane, for the Ising model with h = ±iir/2, 
i.e., fi = —1, on a ladder strip with toroidal boundary conditions. 
Zeros are shown for L x = 150. 



FIG. 4: Complex-temperature phase boundary B u and partition 
function zeros in the u plane, for the Ising model with fj, = —1/2, on 
a ladder strip with toroidal boundary conditions. Zeros are shown 
for L x = 100. 



2u 



1 ± VI + 6u + u 2 



(4.18) 



where the ± signs apply for j = 3, 4, respectively. All of 
these eigenvalues vanish at u = 1, so that Z ~ (u—l) Lx ~ 
(li — l)™/ 2 as u — > 1, i.e., Z has a zero of multiplicity n/2 
at u = 1. 

The boundary is a curve that passes through the 
points u = 0, u = — 1, and l/u = 0, separating the w 
plane into three regions, as shown in Fig. [3] (i) R\, con- 
taining the real interval u > 0, where At£,3 is dominant; 
(ii) i?2, including the real interval u < — 1, where A to is 
dominant; and (iii) i?3, the interior of the loop, including 
the real interval — 1 < u < 0, where Xte,2 is dominant. 
There is an isolated point u = 1 where all four of the 
At^j's, j = 1, ...,4 vanish, and the partition function it- 
self vanishes. The invariance of the locus B u under the 
inversion map u — > l/u is evident in Fig. 03 The inner 
loop is the solution of the equation \Xte,2\ = (At^^l, while 
the outer curve extending to u = ±ioo is the solution of 
the equation |A^i| = |A« i3 |. 



the right. In Fig. [4] we show a plot of zeros for a typical 
value in this range, fi = —1/2. For this case, the arc 
cndpoints in the L x — > oo limit occur at approximately 
u ~ -0.431214±0.3218815z and the real line segment oc- 
cupies the interval 1.051945 < u < 2. These are certain 
zeros of the polynomial 

(u - 2)(64u 7 + 256u 6 - 796u 5 + 



-272u 4 + 352u 3 



135u- 18) 



(4.19) 



that occurs in a square root in the solution of eq. (|4.6[) 
for this case. The line segment that we find on the real 
u axis for —1 < /x < is the analogue of the two line 
segments on the real u axis that we found for this range 
of /i for the model on the square lattice in [l(| (as shown 
in Fig. 6 of that reference). 



G. 



H = e 



F. -l<tx<Q 

As fi increases from —1 toward zero through real val- 
ues, the above-mentioned loop on the fi = — 1 locus B u 
breaks, with its two complex-conjugate arcs retracting 
from u = 0. These arcs cross each other at u = — 1 with 
the outer parts continuing to extend upward and down- 
ward to infinity in the u plane, passing through l/u = 0. 
The boundary B u separates the u plane into two regions, 
Ri to the right, and i?2 to the left, of these semi-infinite 
arcs. The single zero with multiplicity n/2 that had ex- 
isted at u = 1 for ji = — 1 is replaced by a finite line 
segment in the region /i > 1. As fi moves to the right 
from —1 toward zero, the real line segment also moves to 



Here we analyze the complex-temperature phase dia- 
gram for this strip in the case where h is pure imaginary, 
i.e., fj, = e 10 . We show that the endpoints of the unit- 
circle arc on B M , i.e., the Yang-Lee edge singularities, 
have a corresponding feature in B u , namely an endpoint 
of a real line segment that lies in the interval < u < 1 . 
For /i = e l8 the eigenvalues of the transfer matrix consist 
of A = 1 — u as in eq. (|4. 5[) and the three roots of the 
cubic (|4.6p . The coefficients 02 and a\ in this cubic can 
be expressed conveniently as 



«2 



-(1 +u + 2u 



COS! 



and 



ai 



(1-u) 1 + m(2cos6» + 1) 



(4.20) 



(4.21) 
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FIG. 5: Phase boundary B u in the u plane for the Ising model 
with 11 = i, on the ladder strip with toroidal boundary conditions. 
Zeros are shown for L x = 200. 



We find that for /i on the unit circle, the complex- 
temperature phase boundary always passes through the 
points u = 0, u = — 1, and u = ±ioo (the last corre- 
sponding to the single point 1/u = in the plane of the 
variable 1/u = e 4K ). We now prove these results. To 
show that the point u = — 1 is on B u , we observe that 
for u = —1, the eigenvalue given by eq. (|4. 5[) has the 
value A = 2, and the cubic equation for the other three 
eigenvalues factorizes according to 

(A -2) [A 2 +4 cos 2 ((9/2) A + 4] = , (4.22) 

so that these three other eigenvalues are A = 2 and 

(4.23) 



A 



-cos 2 (cV2)±isin(6>/2) \Jl + cos 2 (6>/2) 



All of these have magnitude 2, which proves that the 
point u = —1 is on B u . Indeed, this calculation shows, 
further, that four curves on B u intersect at u = — 1. To 
prove that the point u = is on B u , we first note that for 
this value of u, the eigenvalue Xte.i — 1 — u has the value 
1. We multiply eq. (|4.6[) by u 2 and then take the limit 
u — > 0, obtaining the equation A — 1 = 0. This proves 
the result since we then have two degenerate dominant 
eigenvalues. The same method enables one to conclude 
that the point 1/u = is on B u . 

For any 0^0 mod tt, the locus B u includes a line 
segment that occupies the interval — 1 < u < and also 
occupies part of the interval [0,1). We denote the right- 
hand end of this line segment as Ute,rhe- This right-hand 
endpoint increases monotonically from to 1 as 9 in- 
creases from to it. 

As an illustration of the complex-temperature phase 
diagram for u on the unit circle, we consider the case 
9 = 7r/2, i.e., [i — i, for which the cubic equation (|4.6|) 
takes the form 



(4.24) 

In Fig. O we show the resultant complex-w phase di- 
agram. The boundary has a multiple point at tt = 
where a real line segment intersects a vertical branch of 
the curve on B u . There are three triple points on B Ul 
namely the one at u = — 1 and a complex-conjugate pair 
in the second and third quadrants. The right-hand of the 
real line segment occurs at 

u u ,rhe - 0.746125 (4.25) 

which is the unique real positive root of the polynomial 

u 8 + 2m 7 - 2m 6 + 26u 5 - 48u 4 + 30m 3 - 2m 2 - 2u - 1 

(4.26) 

that occurs in a square root in the exact solution of the 
cubic. B u also includes two complex-conjugate curves 
that extend upward and downward to ±ioo within the 
first and fourth quadrants, passing through 1/u ~ 0. As 
is evident from Fig. the boundary B u separates the u 
plane into four regions: R\ and i?2, extending infinitely 
far to the right and left, and the two complex-conjugate 
enclosed phases separated by the part of the real line 
segment — 1 < u < 0. Qualitatively similar results hold 
for other values of fj, = e 10 with < 9 < tt. In Table 
U we show the values of u r h e ,ti and the corresponding 
values of fc^T/J, denoted as kBT r h e ,t£/ J , as functions of 
9. These are compared with the values for the Ising model 
on the infinite line (with periodic boundary conditions), 
denoted, respectively, as u r h e ,iD and kBT r h e ,iD / J ■ 



TABLE I: Values of 



and k B T rhe /J = K , as a 



function of 9, with = e , for ID (columns 2,3) and the toroidal 
lattice (tt) strip (columns 4,5). We use the notation T = kgT/J. 
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UrheAD 


Urhe,U 




T r h e ,te 













tt/12 


(2 - a/2 + \/3)/4 ~ 0.017 


0.1918 


0.982 


2.422 


tt/6 


(2 - \/3)/4 ~ 0.0670 


0.3315 


1.480 


3.622 


tt/4 


(2 - a/2)/4~ 0.1464 


0.45245 


2.082 


5.044 


tt/3 


1/4 


0.5612 


2.885 


6.924 


tt/2 


1/2 


0.7461 


5.771 


13.658 


2tt/3 


3/4 


0.8846 


13.904 


32.625 


3tt/4 


(2 + a/2")/4~ 0.85355 


0.9346 


25.261 


59.107 


5tt/6 


(2 + x/3)/4 ~ 0.9330 


0.9707 


57.690 


134.725 


— > 7T 


-» 1 


-> 1 


— > oo 


— > oo 



A 3 - (1 + u)\ 2 + u - 2 {i - u 2 )\ + u- 2 (u - iy 



. 



V. EXACT SOLUTION FOR CYCLIC LADDER 
STRIP 

We have also carried out a similar study of B u for the 
ladder strip with cyclic (i.e., periodic longitudinal and 



9 



free transverse) boundary conditions. As before, we take 
the length L x to be even to maintain the bipartite prop- 
erty of the infinite square lattice. Strips of this type (with 
L x > 4 to avoid degeneration) have the property that all 
of the sites have the same coordination number, 3. They 
are thus not expected to exhibit properties that are as 
similar to those of the square lattice as the toroidal lad- 
der strip (which has the same coordination number as 
the infinite square lattice). Furthermore, since the trans- 
verse boundary conditions are free rather than periodic, 
finite-size effects are larger for this lattice than for the 
toroidal strip, which has no boundaries. Thus, we in- 
clude a discussion of the Ising model on this cyclic ladder- 
strip mainly for comparative purposes, but since the strip 
shares fewer similarities with the infinite square lattice, 
our treatment will be more brief than for the toroidal 
strip. 

Since the cyclic ladder strips have odd coordination 
number, the partition function is a Laurent polynomial 
in /i and in x = e~ 2K , rather than u = x 2 = e~ AK . 
(Here we switch notation from that used in our previous 
papers 0, [l(| EL US], usm g x rather than z for e~ 2K 
in order to avoid confusion with the fugacity z = , 
where // is the chemical potential.) In [llT | we proved 
that for lattices with odd coordination number, the fol- 
lowing equality holds (up to a possible overall factor that 
does not affect the zeros): Z(x,—fi) = Z(—x,/.i). This 
theorem implies that the accumulation set of the zeros 
satisfies 



B(x,-fj,) = B{~x,fi) 



(5.1) 



In particular, for real /i, this symmetry, together with 
the inversion symmetry /i — > l//i, means that it suffices 
to consider just the interval < fi < 1. In the basis (|4.1|) 
the transfer matrix is 



/ e 3K+2h 



-3K 



-3K 



V 



-K 



e~ K 

„3K-2h 



\ 



(5.2) 



where c£ denotes "cyclic ladder" . This transfer matrix 
has the same determinant as for the toroidal ladder. The 
partition function is given by 



Z c e - y^(Arfj ) i 



(5.3) 



where \ c t,j, j — 1, ..4 are the eigenvalues of T C £. the 
reduced free energy is / = (l/2)\n(X c e, max ). The eigen- 
values are A c ^i = x~ x l 2 {\ — x 2 ) and the three roots of 
the cubic polynomial that comprises the rest of the char- 
acteristic polynomial of T c g. The reduced free energy is 
/ = (l/2)ln(A rf 

.max ) • 

For h = the phase boundary B x for the g-state Potts 
model and, in particular, the q — 2 Ising case, was an- 
alyzed in [16J. This boundary consists of curves that 



pass through x = 0, x — ±i, and 1/x = with each of 
these points being a multiple point where two branches 
of B x cross each other. The curves separate the x plane 
into six regions, two of which include the real intervals 
x < and x < 0, and the other four of which include 
the intervals on the imaginary axis — oo < Im(cc) < — 1, 
-1 < Im(x) < 0, < Im(x) < 1, and 1 < Im(a;) < oo. 
There is also an isolated zero of the partition function at 
x = — 1 with multiplicity scaling like n. 

Here we focus on the case of nonzero h. For this case 
the two pairs of complex-conjugate curves connecting x = 
with x — ±i each break, and the endpoints move away 
from the real axis as fi decreases from 1 to 0; at the 
same time, the multiple zero at x — —1 is replaced by a 
line segment. One of the reasons for studying this lattice 
strip is to confirm that the singular locus B x again has 
a line segment in the physical ferromagnetic region, just 
as we found for the ID line and the toroidal strip. We 
do, indeed, confirm this, showing the generality of this 
important result. A general property of B x for the Ising 
model on this cyclic strip is that for /x on the unit circle, 
fi = e* e , this boundary B x passes through x — 0, 1/x = 0, 
and x = ±i. For /i ^ ±1, there are line segments on the 
real axis. The value fi = i is especially simple, since 
the invariance of B under fj, — > 1/fi and the symmetry 
for lattices of odd coordination number, (|5.ip . together 
imply that for fj, = i, B x is invariant under x — ► — x 
and hence can be depicted in the u plane. We find that 
for this case /i = i, the right-hand endpoint of the real 
line segment occurs at x r h e — 0.82942, or equivalently, 
u r he — 0.68794, which is the unique positive root of the 
polynomial 



u 6 + 29u 4 



27u z - 4m - 1 



(5.4) 



which occurs in a square root in the exact solution to the 
cubic equation for the eigenvalues of the transfer matrix. 
Since the coordination number of this cyclic lattice is 
intermediate between the value, 2, for the periodic ID 
line and the value 4 for the toroidal ladder strip, one 
expects that the value of u r he = x 2 he at a given value of 
9 would also lie between those for the ID line and the 
toroidal strip. This is verified; we find (see Table [TJ) the 
respective values u r } ie = 0.5, 0.6879, and 0.7461, for the 
ID line, and cyclic toroidal ladder strips. These values 
increase monotonically as the strip width increases and 
can be seen to approach the value of u r h e — 0.8 that we 
infer for the thermodynamic limit of the square lattice 
from our calculations of partition function zeros, to be 
discussed below. 



VI. RELATIONS BETWEEN COMPLEX-u 
PHASE DIAGRAM FOR THE ISING MODEL IN 
ID AND 2D FOR REAL fj, 

In this section we give a unified comparative discussion 
of how our exact results for B u on quasi-lD strips relate 
to exact results for B u in 2D for [i = ±1 and the case 



10 



of real \x in the interval — 1 < fi < 1 that we studied 
earlier in [10j . We first review some relevant background 
concerning the phase diagram for the two cases where 
this diagram is known exactly for the 2D Ising model, 
namely jj, = 1 (h = 0) and fj, = — 1 (/i = ?7r/2) 

A. n = 1 

The complex-u phase boundary ,B M for the square- 
lattice Ising model is the image in the u plane of the 
circles @ 

|ar±l|=V2, (6.1) 
namely the limagon (Fig. lc of [20() given by 
Re(u) = 1 + 2V2 cos(w) + 2 cos(2w) 

Im(w) = 2V2 sin(w) + 2sin(2w) (6.2) 

with — 7r < oj < 7r. The outer branch of the limagon in- 
tersects the positive real-w axis at upm-afm = 3 + 2^/2 
(for to = 0) and crosses the imaginary-w axis at u = ±(2+ 
V3)i (for lu = ±5tt/12). The inner branch of the limagon 
crosses the positive real axis at upm-fm — 3 - 2V2 (for 
u) = tt) and the imaginary-u axis at u — =f(2 — \/3 )i 
(for u> = ±ll7r/12). The limagon has a multiple point at 
u = — 1 (for u> — ±57r/4) where two branches of B u cross 
each other at right angles. When u = — 1, there are also 
branches of the limagon passing through Im(w) = ±2\/2i 
(for lu = ±7r/2). The boundary Z? u separates the u plane 
into three phases, which are the complex extensions of the 
physical PM, FM, and AFM phases. Since the infinite- 
length strips are quasi-lD, the Ising model has no finite- 
temperature phase transition on these strips, and is crit- 
ical only at T = 0. Thus, for h = 0, the boundary B u 
passes through u = and 1/u = 0. However, just as 
for the square lattice, for the toroidal strip the boundary 
B u separates the u plane into three regions, as is evi- 
dent in Fig. [1] One can envision a formal operation on 
the boundary curve B u for the toroidal ladder strip that 
transforms it into the B u for the 2D lattice, namely to 
move the crossing at u = to upM-FMi which, owing to 
the u <-> 1/u inversion symmetry, automatically means 
that the part of the boundary B u that goes to infinity in 
the u plane is pulled back and crosses the real axis at the 
inverse of this point, viz., upm-afm ■ 

The complex-w phase boundaries B u of the Ising model 
on both the infinite-length ID line and on the infinite- 
length ladder strip with toroidal or cyclic boundary con- 
ditions have the property that they pass through u = — 1 
and, for the toroidal and cyclic ladder strips this is again 
a multiple point on B u , just as it is in 2D. For the toroidal 
strip, the point u = —1 is an osculation point, where two 
branches on B u intersect with the same tangent, whereas 
for the square lattice the branches cross at right an- 
gles. Other similarities include the fact that, e.g., for 
the toroidal strip, B u crosses the imaginary u axis at two 



pairs of complex conjugate points that are inverses of 
each other, namely u — ±(a/2 + l)i and u = ±(\/2 — 
These points are in 1-1 correspondence with the points 
u = ±(2 ± i/3 )i where B u crosses the imaginary- it axis 
for the square lattice. 

B. n = -l 

The phase boundary for the Ising model with fj, = — 1 
on the square lattice was determined in [9( and consists 
of the union of the unit circle and a line segment on the 
negative real axis: 

B„( M = -1): {|u| = 1} 

U {-(3 + 2V2) < u < -(3- 2V2)} . (6.3) 

It is interesting that the endpoints of this line segment 
are minus the values of upm~fm and upm-afm = 
l/upM-FM on the square lattice. The point u = — 1 
is a multiple point on B u where the unit circle \u\ = 1 
crosses the real line segment at right angles. The latter 
feature is matched by the locus B u for the toroidal ladder 
strip, as is evident in Fig. [3l 

C. < u < 1 

In [l(| it was found that as h is increases from 0, i.e., 
as n decreases from 1 to 0, the inner loop of the limagon 
immediately breaks open at u = upm-fm, forming a 
complex-conjugate pair of pr ong endpoints u e , it* that 
retract from the real axis. In hfj we used calculations of 
complex-M partition function zeros together with analy- 
ses of low-temperature, high-field series to determine the 
locations of these arc endpoints and the values of the ex- 
ponents a' e , (3' e , and 7^ describing the singular behavior 
of the specific heat, magnetization, and susceptibility at 
these prong endpoints. 

In contrast, the PM-AFM critical point does not dis- 
appear. For the antiferromagnetic sign of the spin-spin 
coupling, J < 0, as H increases, the Neel temperature T/v 
decreases, or equivalently, — K c = \ J\/(kpTN) increases, 
and hence also uzm-afm increases from its value of 
3 + 2V2 at H = (where the notation ZM follows [13]). 
As H increases sufficiently, there is a tricritical point, 
and when it increases further to — cJ = c\J\, where c 
denotes the coordination number, the Neel temperature 
is reduced to zero. This means that /3zm-afm —> °o- 
Thus, asymptotically as h — > 00, —KzM-AFhi/h — > c, 
i.e., the right-hand side of the boundary B u moves out- 
ward to infinity like u ~ /i -1 / 2 as /1 — * 0. 

With the replacement of the finite-temperature PM- 
FM critical point u = u c by the zero-temperature crit- 
ical point u — 0, the boundary B u for the Ising model 
on the infinite-length limit of the toroidal ladder strip 
reproduces this feature of the model in 2D, viz., imme- 
diate breaking of the loop, as is evident in Fig. [21 The 
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corresponding boundary B z for the infinite-length limit 
of the cyclic strip also immediately breaks apart from the 
zero-temperature PM-FM critical point at u = 0. 



D. -1 < (j, < 

For /i in the interval — 1 < (i < 0, in the 2D case, 
the partition function zeros calculated in [1(3] exhibited 
patterns from which one could infer that in the thermo- 
dynamic limit the resultant accumulation locus B u ex- 
hibited curves and two line segments, one on the posi- 
tive, and one on the negative real u axes. (There were 
also some zeros that exhibited sufficient scatter that one 
could not make a plausible inference about the asymp- 
totic locus in the thermodynamic limit.) One may thus 
ask if we obtain qualitatively similar behavior with these 
exact closed-form solutions for the model on quasi-one- 
dimensional strips. For both the ID line and the toroidal 
ladder strip with this range of fj,, we find a line segment on 
the positive real axis (cf. eq. (|3.8I and Fig. HJ, in agree- 
ment with this feature that we had obtained in 2D. As 
could be expected, the quasi-lD strips do not reproduce 
all of the features that we found for 2D. For example, 
neither the ID line nor the toroidal strip exhibits a real 
line segment on the negative real axis, either for the case 
fj, = — 1 or the range — 1 < fj, < 0, where we did find such 
a line segment in 2D. 



VII. COMPLEX-u PHASE DIAGRAM AND 
ZEROS OF THE PARTITION FUNCTION FOR 
THE SQUARE LATTICE WITH fi = e iB 

A. Motivation and Exact Results for 8 = and 

8 = 7T 

In this section we present our calculations of complex-w 
zeros of the partition function of the square-lattice Ising 
model for imaginary h, i.e., fi = e~ 2h = e l9 with < 
9 < it. Owing to the invariance of the model under the 
inversion (|2.2[) . it suffices to consider this half-circle. This 
study is a continuation of our earlier investigation in [t| 
Hoj of the complex-w phase diagram of the model for real 
nonzero external magnetic fields (hence < (i < oo) 
and the subset of complex h of the form h = h r ± iir/2 
yielding negative real /k, and thus covering the interval 
— oo < fj, < 0. 

As noted above, in studying the complex-u phase dia- 
gram, it is natural to consider paths in the fj, plane that 
connect the two values for which this phase diagram is 
exactly known, namely /i — 1 (h = 0) and fi = — 1 
(h = in/ 2). In [l(| we considered the path defined by 
the real interval —1<(1<1. Here we concentrate on 
the other natural path, namely an arc along the unit cir- 
cle \i = e l9 with < 9 < n. For \i on this unit circle we 
have mentioned above that the boundary B u is invariant 
under complex conjugation. One motivation for study- 



ing the complex-u zeros of the partition function for fi on 
the unit circle is that the latter locus is precisely where 
the zeros of the complex-/! zeros of the partition function 
occur for physical temperatures in the case of ferromag- 
netic couplings. Hence, our results in this section consti- 
tute an investigation of the pre-image in the u plane, for 
the square-lattice Ising model, of points on the Yang-Lee 
circle. Indeed, just as we found with exact results on 
infinite-length quasi-lD strips, our calculations of parti- 
tion function zeros for the model in 2D will lead us to 
the inference that in the thermodynamic limit the locus 
B u for /i = e 16 with 9^0 mod tt contains a line segment 
extending into the physical ferromagnetic region with a 
right-hand endpoint u r i ie that corresponds precisely to 
the temperature for which the points fi = e ±ie are the 
endpoints (Yang-Lee edges) of the arc of the unit circle 
= 1 comprising B^. A convenient feature for the study 
of the complex-u phase diagram for the square-lattice 
Ising model with n = e 10 is that the phase boundary B u 
remains compact throughout the entire range of 9. This is 
in contrast to the situation for the real path — 1 < /a < 1. 
In that case, as was discussed in [loj , the phase boundary 
separating the phase where the staggered magnetization 
M st vanishes identically from the AFM phase where M st 
is nonzero moves outward to complex infinity as fi — » 
and then comes inward again as fi passes through and 
approaches (J, = — 1. (It should be noted that although 
B u is compact for n = e 16 with 9 £ K for the square lat- 
tice, this is not the case with the triangular lattice. On 
that lattice, for both of the exactly solved cases 9 = and 
9 = tt, the locus B u contains the respective semi- infinite 
line segments — oo < u < —1/3 and — oo < u < —1/2 

In our previous work ,Tc|, we tested several different 
types of boundary conditions including doubly periodic 
(toroidal, TBC) and helical boundary conditions (HBC). 
The latter are periodic in one direction, say L x , and heli- 
cal in the other, say L y . We found that helical boundary 
conditions yielded zeros that showed somewhat less scat- 
ter for general fi and were closer to the exactly known loci 
B u for the cases // = ±1 than the zeros obtained with pe- 
riodic boundary conditions. This can be interpreted as a 
consequence of the fact that for toroidal boundary con- 
ditions, the global circuits around the lattice have length 
L x and L y , while for helical boundary conditions, while 
the circuit in the x direction is still of length L x the 
one in the y direction is made much longer, essentially 
L x L y . For the present work we have again made use of 
helical boundary conditions and also a set of boundary 
conditions that have the effect of yielding zeros that lie 
exactly on the asymptotic loci B u for the exactly known 
cases /j, — ±1. These are defined as follows. We con- 
sider two L x x L y lattices, with the x direction being 
the longitudinal (horizontal) and and the y direction the 
transverse (vertical) one. We impose periodic longitudi- 
nal boundary conditions and fixed transverse boundary 
conditions Specifically, we fix all of the spins on the top 
row to be + while those on the bottom row alternate in 
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sign as (H 1 — ...). For the second lattice, we impose 

spins on the top and bottom rows that are minus those 
of the first lattice; that is, all spins on the top row are — , 

while those on the bottom are ( — I h ...)■ Together, 

these yield a partition function that is invariant under 
the h — > — h symmetry. We denote these as symmetrized 
fixed boundary conditions (SFBC). In passing, we note 
that if one used only the first lattice, the corresponding 
boundary conditions would correspond to set A of [26| . 
For h = this set was shown to yield zeros that lie exactly 
on the circles (16.11) . For our present work, the boundary 
conditions of [26| would not be appropriate, since they vi- 
olate the h — * — h symmetry and hence also the /i — > 1/fi 
symmetry of the infinite square lattice. In turn, this vi- 
olation would have the undesirable consequence that for 
/i = e %e , the set of zeros would not be invariant under 
complex conjugation and zeros that should be exactly on 
the real-w axis would not be. We have found that the 
symmetrized fixed boundary conditions yield zeros with 
somewhat less scatter than helical boundary conditions, 
and therefore we concentrate on the former in presenting 
our results here. 

For the analytic calculation of the partition function, 
we again use a transfer matrix method similar to that 
employed in our earlier paper (Io| . In that work we per- 
formed a number of internal checks to confirm the ac- 
curacy of the numerical calculations of the positions of 
the zeros of the partition function. Since for our present 
study we are performing calculations of partition func- 
tions and zeros for considerably larger lattices than we 
used in [icj |. we have paid special attention to guaran- 
teeing the accuracy of the numerical solution for these. 
Among other things, we now use the rootsolver program 
called MPSolve 27] to augment the internal rootsolvers 
in Maple and Mathematica. 

We show our results for the complex-it zeros of the 
Ising model partition function on L x x L y sections of the 
square lattice with /i = e l9 in Fig. [6] for various values 
of 9 in the range < 9 < tt. These zeros were calculated 
with the symmetrized fixed boundary conditions defined 
above. The curve (|6.2|) for 9 = and the curve and line 
segment (|6.3|) for 8 = tt represent exact results. A more 
detailed view of the inner region near u = is shown in 
Fig. [71 We present a detailed view of the zeros in the 
inner central region for tt/2 < 8 < tt in Fig. [5) The 
zeros presented for 8 ^ tt/2 were calculated on 12 x 13 
lattices. As discussed further below, we devoted a more 
intensive study to the value 8 = 7r/2, i.e., /i = i, and for 
this case we calculated the partition function and zeros 
for L x x L y lattices with sizes L x and L y ranging from 12 
to 16. We show the results for this /j, = i case separately 
in Fig. [9] Concerning exact results, for visual clarity, 
in Figs. [7J and [5] for the case 8 = tt, we show only the 
right-hand endpoint of the real line segment (|6.3[) on B u 
at u = — (3 — 2^2) (indicated by the symbols <3 and 0, 
respectively). 

As 8 increases from zero, we observe a number of in- 
teresting features of the complex-u zeros of the partition 
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FIG. 6: Complex- ?i zeros of the Ising model partition function on 
sections of the square lattice for = e with 8 = 7r/4, tt/2, 3n/4, 
and ll7r/12. See text for details of calculation. The exact phase 
boundaries for 6 = and 8 = tt are also shown. 
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FIG. 7: A closer view of the complex-tt partition function zeros in 
the region near u = for several values of 8 in the range < 8 < n. 
For reference the figure also shows, as exact results, the inner part 
of the curve forming the limagon i l(i.2[ i for 8 = and the unit circle 
on B u for 8 = tt. For visual clarity, for 8 = tt, we show only the 
right-hand endpoint, u = —(3 — 2v / 2), of the real line segment on 
B u (indicated with the symbol <). The lattice size is 12 x 13 except 
for 8 = it/2, for which we show results with L x and L y up to 16. 



function. One important general feature is that for 9 =/= 
mod tt, zeros occur on the real axis, extending over an 
interval from the point where the inner loop of B u is in- 
ferred to cross this axis, to a right-hand endpoint u r h e 
that increases as 9 increases in the interval < 9 < tt. 
We infer that in the thermodynamic limit (i) these ze- 
ros merge to form a real line segment on B u and (ii) the 
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FIG. 8: A closer view of the complex- u partition function zeros 
in the inner central region for a several values of 6 in the range 
7r/2 < 8 < 7T. Lattice size is 12 X 13, as in Fig. \E\ 
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FIG. 9: Complex-u zeros of the Ising model partition function 
with = e 1 "^ 2 = i on L x X L y sections of the square lattice with 
aspect ratio ~ 1 and L x and L y varying from 12 to 16. 

right-hand endpoint 



corresponds to the temperature T r /, e at which the circular 
arc comprising B^ has endpoints at /i — e ±l6 . For infi- 
nite temperature, K r h e — 0, this endpoint (the Yang-Lee 
edge) occurs at 9 = tt and as the temperature decreases, 
the endpoints of the circular arc on 2? M moves around to 
progressively smaller values of 9. As T decreases to the 
critical temperature T c = Tpm-fm for the onset of fer- 
romagnetic long-range order, 6 1 — > 0, closes to form 
the unit circle = 1, and for lower temperatures it 
remains closed. The property that u r h e and the corre- 
sponding temperature T r h e increase monotonically with 



9 in the range < 9 < tt is equivalent to the property 
that the complex-conjugate endpoints of the circular arc 
comprising £? u (i.e., the Yang-Lee edge) at 9 increases 
monotonically from 9 = at T = T c to 9 — > tt as T — > oo. 
This thus establishes a 1-1 correspondence between the 
right-hand endpoint u r he of the real line segment for a 
given 9 and the temperature at which this 9 is the value 
of the endpoint of the circular arc on B^. The limit 9^0 
involves special behavior, in that this real line segment 
shrinks to zero and disappears. The limit 9 — > tt is also 
special; again, the line segment on the positive real u axis 
disappears in this limit and is replaced by a single zero at 
u = 1 with multiplicity n/2, where n denotes the number 
of sites on the lattice. This zero at u = 1 gives rise to the 
term (1/2) ln(u — 1) in the reduced free energy at /i = —1 

aai. 

Turning on a finite (uniform) magnetic field, whether 
real or complex, does not remove the PM-AFM phase 
transition that occurs for sufficiently large negative K = 
f3J. It follows that the outer loop on B u cannot break. 
As discussed in [lOj, this can be shown via a proof by 
contradiction. Assume that this outer loop on B u did 
break; then one could analytically continue from the re- 
gion around K = 0, i.e., u — 1, where the staggered 
magnetization M st vanishes identically, to the physical 
AFM phase where M st is nonzero, and similarly to the 
complex-it extension of this AFM phase, which would 
be a contradiction. We find that as 9 increases from 
to 7r, the zeros that form the outer loop of B u in 
the Re(it) > half-plane move monotonically inward 
toward the unit circle |u| = 1, which they form for 
9 — > tt. We infer that in the thermodynamic limit, (i) 
the right-most crossing on B u decreases monotonically 
from u = upM^AFAi = (3 + 2^/2) ~ 5.83 to u = 1 as 
9 increases from to tt; and (ii) the upper and lower 
points where the outer loop of B u crosses the imaginary 
u axis move monotonically inward from u ~ ±(2 + a/3)« 
to u = ±i. From inspection of the actual zeros that we 
calculate for various values of 9, we infer the following 
approximate maximal values of u at which B u crosses 
the positive real u axis: u ~ 5.5 for 9 — tt/A, u ~ 4.6 for 
9 = 7r/2, and u ~ 3.3 for 9 = 3tt/4. 

Recall that for nonzero real h, a theorem [2|| guaran- 
tees that the free energy for the ferromagnet is analytic 
for all temperatures, which means that B u must break 
and retract from the real axis in the vicinity of what was, 
for h = the PM-FM phase transition point, upa/-fa/ 
[25| . In contrast, the results of 0, [2{| allow the free en- 
ergy to be non-analytic as H is varied at constant (3 or 
/3 is varied at constant H (i.e., in both cases, as h is var- 
ied) if Re(h) — 0, that is, h is pure imaginary, which is 
the situation that we consider here. Our most detailed 
study of the partition function zeros, for fj, = i (see Fig. 
[5]) is consistent with the conclusion that the inner loop 
on B u does not break but remains closed. This conclu- 
sion is also consistent with our results for other values of 
9. In making this statement, we note that the fact that 
the zeros on the right-hand side of the inner loop, calcu- 
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lated on finite lattices, do not extend all the way in to 
the real axis, does not constitute evidence of a break in 
this loop in the thermodynamic limit. For example, even 
for the exactly solved case fj, = 1, the zeros calculated on 
finite lattices also do not extend all of the way down to 
the real axis. In this context, we also remark on our ex- 
act results for quasi-lD strips; on both the toroidal and 
cyclic ladder strips, for /i = e 1 ® , as 9 is increased from 
to 7T, the loop on B u that passes through the critical 
point (at u = 0) remains intact and unbroken. (Note 
that the point at which this loop crosses the real axis 
for these quasi-lD infinite-length ladder strips remains 
at u = as 9 increases from to 7r, while for the model 
in 2D the inferred crossing point of the inner loop moves 
gradually to the left as 9 increases through this range for 
the square lattice.) 

The details of the pattern of zeros in the complex-it 
region that includes the real interval — 1 < u < are 
complicated, and there is significant scatter of some of 
these zeros. Consequently, we do not try to make further 
inferences about the form of the complex-u boundary B u 
in this region in the thermodynamic limit. As an ex- 
ample of the kind of feature that might be present in 
this limit, one can discern some indication of possible 
triple points at it ~ — 0.7± 0.2i and u ~ —0.9 ± 0.5i. A 
complex-conjugate pair of triple points is, indeed, present 
in our exact solution for B u on the toroidal ladder strip 
with fi = i, as shown in Fig. The scatter of ze- 
ros in this region raises the question of whether some 
part of B u might actually fill out 2-dimensional areas 
rather than being one-dimensional (comprised of curves 
and possible line segments) in the thermodynamic limit. 
For the 2D Ising model in zero field it is easy to see that 
complex-it zeros generically fill out areas if the spin-spin 
exchange constants in the x and y directions are unequal, 
but this is not directly relevant here, since we only con- 
sider the model with isotropic couplings. For isotropic 
couplings, this area behavior happens for a heterop olyg - 
onal Archimedean lattice, namely the 4-8-8 lattice [28j |. 
and here again, the origin of this is obvious from the ex- 
act form of the free energy (see eq. (6.5) and Fig. 7 of 
[111 ) . One can fit curves or line segments to many of the 
zeros in Fig. [6] As for the region where the zeros show 
scatter, our results are not conclusive, and we do not try 
to make any inference about whether or not some set of 
these zeros might merge to form areas in the thermody- 
namic limit. 

Several other aspects of the real zeros are of interest. 
First, we find that as 9 increases from 0, there are real 
zeros not just to the right of the extrapolated point where 
the inner loop on B u crosses the real axis, but also to the 
left of this point. Indeed, we find that for < 9 < tt, 
there are zeros on the negative real axis. As 9 — > tt, these 
occur in the interval of eq. dO}, i.e., -(3 + 2\/2) < u < 
— (3 — 2^2). On the finite lattices that we have studied, 
we also have found complex-conjugate pairs of zeros that 
are close to the zeros on the real axis. 

The value 9 = 7r/2, i.e., /x = i, is the middle of the 



range under consideration here, and we have devoted a 
particularly intensive study to it. In addition to the gen- 
eral plots comparing the zeros for this value of 9 with 
those for other values of 9, we show the zeros for this case 
alone in Fig. [9l The zeros shown in this figure were cal- 
culated for several different lattice sizes with aspect ratio 
L y /L x ~ 1 and L x ranging from 12 to 16. We display the 
zeros for different lattice sizes together to see lattice size- 
dependent effects. As is evident from the figure, much 
of the locus comprised by these zeros is largely indepen- 
dent of lattice size for sizes this great. These calculations 
illustrate our general description of the properties of B u 
above. From these zeros we infer that for fj, = i, in the 
thermodynamic limit, (i) the complex- it phase boundary 
B u crosses the real axis at u ~ 4.6 and the imaginary 
axis at 11 ~ ±2.7i; (ii) there is an inner loop on B u which 
is likely to pass through u = 0, although there is some 
decrease in the density of complex zeros in the vicinity 
of this point; (iii) the locus B u exhibits a line segment 
on the real axis that extends from a right-hand endpoint 
u r h e — 0.8 leftward with components along the negative 
real axis; and although the details of this line segment 
at intermediate points cannot be inferred with certitude, 
the left-most endpoint occurs at uihe — —4.5; (v) the 
phase boundary B u thus appears to separate the it plane 
into at least four regions: (a) the AFM phase and its 
complex-it extension, which occupies values of u going 
outward to complex infinity, (b) the interior of the outer 
loop; (c) and the complex-conjugate pair of regions in- 
side of the inner loop, which seems to be divided into an 
upper and lower part by the real line segment inside this 
loop. As regards item (i), the points at which the outer 
loop of B u crosses the imaginary it axis are consistent, to 
within the accuracy of our calculations, with being equal 
to ±(1 + y/S)i. 



B. Connections with Results on Quasi-lD Strips 

for ^ = e 10 

With appropriate changes to take account of the 
change in dimensionality, we can relate these features 
to our exact results on quasi- ID strips. For these strips 
we found that the locus B u includes a line segment on 
the positive real axis in the physical ferromagnetic inter- 
val as 9 increases above zero. In the ID case, we found 
the simple result u r h e = sin 2 (6*/2) in eq. (|3.9p for the 
right-hand endpoint of this line segment. For the toroidal 
and cyclic ladder strips we illustrated, e.g. for 9 = it/2 
(fi = i), how it is determined as the root of the respec- 
tive polynomials, eqs. (I4.26P and (|5.4p . that occur in the 
solution of a relevant cubic equation for the eigenvalue of 
the transfer matrix. In Table [T] we showed the values of 
Urhe for the ID line and the toroidal strip, together with 
the corresponding values of the temperature T r h e , as a 
function of 9 = 9 e . This table shows how u r he and T r h e 
increase as 9 increases above and approaches it. We 
also noted how, for a given value of 9, T r } le increases as 
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one increases the strip width. This is physically under- 
standable, since a given value of the angle 9 corresponds 
to a higher temperature and hence larger u r h e as the strip 
width increases because that increase fosters short-range 
ferromagnetic ordering. 

Another property of the zeros that can be related to 
our exact results on quasi- ID strips is the part of the 
line segment extending to the left of the point where the 
inner loop appears to cross the real axis and, indeed, ex- 
tending to negative real values. For the quasi-lD strips, 
these intervals are the same, since the ferromagnetic crit- 
ical point is at u — 0. For the ID line case, the locus B u , 
which is u < sin 2 (0/2), includes the semi-infinite line seg- 
ment u < 0. For the toroidal ladder strip, we find that 
for any 9 in the interval < 9 < 7T, B u includes the seg- 
ment — 1 < u < as well as the portion on the positive 
real axis discussed above. So there are again similari- 
ties with respect to this feature as regards the results 
for the strips and for our zeros calculated on patches of 
the square lattice. In future work it would be of interest 
to calculate complex-u zeros of the Ising model partition 
function with imaginary h on d = 3 lattices, as was done 
for real h in Ref. [29fl . 



VIII. CONNECTION OF SINGULAR 
BEHAVIOR OF THE ZERO DENSITY IN THE n 
AND u PLANES 



For the ferromagnetic Ising model, studies have been 
carried out of the singularity at the endpoint of the circu- 
lar arc Bfj, as 9 — > 9 e (Yang-Lee edge) and the associated 
density of zeros g(9) in the original papers [4| and in 
works including those by Griffiths, Fisher and collabo- 
rators, and Cardy [3(|-[33j]. Kim has suggested that the 
singularity at an edge of a locus of zeros in the u plane 
is equivalent to the Yang-Lee edge singularity jjyj. For 
the case of d = 2 dimensions, we can show this equiv- 
alence using conformal field theory methods. We recall 
that for a conformal field theory indexed by (relatively 
prime) positive integers p and p' , the central charge c is 



c= 1-6 
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with scaling dimensions 



/\2l 



[(pr-p's)' 2 - (p-p') 

hr - s = w ' 

where 1 < r < p — 1 and 1 < s < p' — 1 . The requirement 
of a single scaling field (other than the identity) leads 
uniquely to the identification of the conformal field the- 
ory as which is non-unitary, with central charge 
c = —22/5 [33]. The scaling dimension for the single 
non-identity field is r\ — 4/i 12 = —4/5. From this and 
the standard relation rj = d + 2 — 2yf l , where y^ is the 
magnetic exponent, it follows that y^ = 12/5. Since the 



theory has only one relevant scaling field, the thermal 
exponent is the same, i.e. 



yt=Vh = y ■ ( 8 - 3 ) 

Substituting the value of rj into the sca ling relation (|3.24[) 
with d = 2 yields the result a = —1/6 [32t l33|. The equal 
thermal and magnetic exponents in eq. (|8.3p determine 
all of the rest of the exponents for this critical point, 
which include 
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(These results have also been independently and simul- 
taneously obtained in this manner by B. McCoy.) The 
fact that the conformal field theory has only a single non- 
identity operator and equal thermal and magnetic expo- 
nents leads to the conclusion that the exponent 1 — a' e 
describing the singular behavior of the density of zeros 
at the right-hand endpoint u e = u r he of the positive real 
line segment on B u corresponding to a value of 9 with 
p = e 10 , < 9 < 7T, is the same as the exponent a de- 
scribing the singular behavior of the density of zeros at 
the complex-conjugate endpoints of the circular arcs on 

A parenthetical remark may be of interest here. We 
recall that a liquid-gas phase transition may be mod- 
elled as a lattice gas, and the latter may, in turn, be 
mapped onto a ferromagnetic Ising model. In this con- 
text the Yang-Lee circle theorem states that the zeros of 
the partition function occur on the unit circle \z\ = 1 in 
the complex plane of the fugacity, z — e' 3 ' 1 (where fj,' is 
the chemical potential) and, in the thermodynamic limit, 
form an arc of this unit circle extending from 9 = tt over 
to complex-conjugate endpoints at z — e ±ie " [4]. For 
closed-form approximations to the equation of state of a 
liquid-gas system such as that of van der Waals, the den- 
sity of Yang-Lee zeros has been calculated [3^, H3| ■ For 
intermolecular potentials with a repulsive hard core, the 
expansion of the reduced pressure pv^j '(fcgT) in powers 
of fugacity, pvo/ftsT) = Yl'jLtbezr, exhibits alternat- 
ing signs, indicating a singularity on the negative real 
z axis [36|. This latter singularity has been shown to 
be equivalent to the Yang-Lee edge singularity [H, Hi| . 
This is a somewhat different equivalence than the one 
discussed here, since it relates the Yang-Lee singularities 
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at z = e ±ie " to one on the negative real z axis, whereas 
the relation discussed here is between the Yang-Lee sin- 
gularities and a singularity on the (positive) real u axis. 

Now consider a switch from the imaginary values of h 
relevant for the Yang-Lee edge singularity to real values 
of h. These lead to the complex-conjugate arcs on B u 
with arc (prong) endpoints u e and it* that retract from 
the position of what was the critical point at u = u c 
(for h = 0) as \h\ increases from zero Again, the 

fact that the conformal field theory has only a single 
non-identity operator and equal thermal and magnetic 
exponents leads us to the further conclusion that the ex- 
ponent 1— a' e describing the singular behavior of the zeros 
at the ends of the complex-conjugate arcs in the complex 
u plane at u e and u* (cf. eq. (|3.15|) with u s — u e ) is 
the same as the exponent a. With a' e = 7/6, this implies 
that for the 2D Ising model, g(u) thus diverges at these 
arc endpoints with the exponent 

l-a' e = a = ~ (2D), (8.8) 

We also conclude that the exact values of the exponents 
a' e , /3 e , and "f' e for the specific heat, magnetization, and 
susceptibility in eqs. (|8.5[) - (|8.7p apply to the arc end- 
points at u e and «*. These exact values agree very well 
with the numerical values that we obtained in [l(| from 
our analysis of low-temperature, high-field (i.e., small- 
it, small- jtt) series (see Table I of [l(|). These values 
had been suggested in [34| based on the assumption that 
lit = Uh = 12/5 for the endpoint of a locus of zeros in 
the u plane. Here we have proved the equivalence using 
conformal field theory methods. 

In [10l | we also studied complex values of h correspond- 
ing to negative \i in the real interval — 1 < fx < 0. For the 
solvable case /i = — 1 one knows the exponents a' s and 
f} s exactly at various singular points, and in [9], from 
analyses of series, we obtained the exponent 7^ and in- 
ferred exact values for this exponent also. These singu- 
lar points at \i = — 1 include the multiple point u = — 1, 
the left- and right-hand endpoints of the line segment 
uthe = -(3 + 2^2) and u rhe = l/u£ he = -j.3-2i/2), and 
also the point u = 1. For reference, in [9J] we obtained 



a' e = 1, /? e = -1/8, and i e = 5/4 at u e = -(3 - 2^2), 
a' s = (Ch finite), (3 S = 1/2, and 7^ = 1 at u = —1, and 
a' s = (C H finite), & = -1/4, and j' s = 5/2 at u = 1 
(see Table 4 of that paper) , where the results for a' and 
[3 were exact and the results for 7' were inferred from 
our analysis of series. (Exponents at u r he — l/u£h e are 
related by the u — > l/u symmetry.) For /i close, but not 
equal, to /i = —1, the line segment on the negative real 
axis shifted slightly, and there appeared a new line seg- 
ment on the positive real axis extending inward from the 
right-most portion of the boundary B u . We also studied 
these singular exponents via series analyses in [io| . In 
future work it would be worthwhile to understand better 
the values of the exponents describing these singularities 
for negative real /i. 

IX. CONCLUSIONS 

In this paper we have studied properties of the Ising 
model in the complex u plane for nonzero magnetic field. 
We used exact results for infinite-length quasi-lD strips 
to provide insights into the previous results that we had 
obtained in [9j,ll0|. We also studied the case of complex h, 
in particular, the case of imaginary h, for which \i = e 10 . 
We used both exact results on strips and partition func- 
tion zeros to analyze the phase diagram in the u plane 
for this range of /1. One important result that we found 
was that the boundary B u contains a real line segment 
extending through part of the physical ferromagnetic in- 
terval < U < 1, with a right-hand endpoint u r h e at 
the temperature for which the Yang-Lee edge singularity 
occurs at fj, — e ±l6 . We also used conformal field the- 
ory arguments to relate the singularities at u r h e and the 
Yang-Lee edge. 
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